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The scattering of bright quantum solitons at barrier potentials in one-dimensional geometries is 
investigated. Such protocols have been predicted to lead to the creation of nonlocal quantum super- 
positions. The centre-of-mass motion of these bright matter-wave solitons generated from attractive 
Bose-Einstein condensates can be analysed with the effective potential approach. An application 
to the case of two particles being scattered at a delta potential allows analytical calculations not 
possible for higher particle numbers as well as a comparison with numerical results. Both for the 
dimer and a soliton with particle numbers on the order of N = 100, we investigate the signatures 
of the coherent superposition states in an interferometric setup and argue that experimentally an 
interference pattern would be particularly well observable in the centre-of-mass density. Quantum 
superposition states of ultra-cold atoms are interesting as input states for matter-wave interferom- 
etry as they could improve signal-to-noise ratios. 



I. INTRODUCTION 

The experimental realisation of macroscopic quantum 
superpositions is a challenge of current research. Be- 
sides opening up a new realm for the study of funda- 
mental quantum effects there would be possible appli- 
cations in interferometry and quantum information sci- 
ence. As quantum superposition states are very sensitive 
to decoherence the creation of these fragile objects is a 
difficult task. So far experimental achievements include 
the realisation of mesoscopic superposition states on the 
single-atom level [JJ , a six-atom Schrodinger-cat state [5] 
and superposition states in radiation fields [3]. Further 
systems of interest for the creation of many-particle en- 
tanglement are cavity quantum optomechanical systems 
[3] and mesoscopic superpositions of internal degrees of 
freedom. In the field of ultra-cold atoms mesoscopic su- 
perposition states are important for matter-wave inter- 
ferometry as they allow improved signal-to-noise ratios 
[5]. Many-particle entanglement in macroscopic systems 
enables investigations of the quantum-classical boundary 
and decoherence mechanisms [§]. Therefore larger in- 
terferometric objects are required. While interferences 
have been observed with fullerenes pj [5] and large or- 
ganic molecules [5] the creation of a macroscopic non- 
local Schrodinger-cat state of ultra-cold atoms still is 
a challenge. Recent proposals have been the creation 
of a nonlocal state of a Bose-Einstein condensate sus- 
pended between the two wells of a two-site optical lat- 
tice [10] or via scattering of bright matter-wave solitons 
(BMWS) [UJIU. Single BMWS [TJ and soliton trains 
[HI IT5] have already been realised experimentally in low- 
dimensional attractive Bose gases. Theoretical investi- 
gations of one-dimensional attractive Bose-gases [16j Hi] 
comprise the investigation of BMWS collisions [JH] and 
collision-induced entanglement [19] , 
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In one-dimensional geometries the scattering of a 
BMWS or a tightly bound dimer off a barrier poten- 
tial like a laser focus has been predicted to lead to 
the creation of a nonlocal quantum superposition of all 
the atoms being located to the left and to the right of 
the scattering potential. In a Fock space representation 
with states I^l^r), where (n^) denotes the particles 
left (right) of the barrier, this would correspond to the 
NOON-state [5D] 

|*NOON> = -^(|iV,0} + e^0,iV)) (1) 

where (f> is the relative phase between the two parts of 
the superposition. Experimental parameters for this pro- 
posed realization of the ideal of a nonlocal quantum su- 
perposition have been estimated in [TTj : for a BMWS 
consisting of N = 100 7 Li atoms with a typical size of 
about 1 /im a spatial separation of about 100 /im is found 
to be reachable. The initial wave-function should be pre- 
pared carefully such that the beam-splitter acts on the 
soliton as a whole and not on each single atom. Using 
the effective potential approach from [11] it is possible 
to investigate the dynamics of the centre of mass (CoM). 
While in Ref [11] the effective potential approach has 
been applied to broader potentials, here it will be applied 
for smaller potentials which is more interesting from the 
point of view of a theoretical physicist: for low particle 
numbers, it allows to perform Rayleigh-Schrodinger per- 
turbation theory beyond the leading order of the effective 
potential investigated in Ref. [TT]. This thus renders a 
detailed test of the effective-potential approach possible. 

Two-particle bound states are interesting both exper- 
imentally [3TJ and theoretically p2!?M25| . An application 
of the effective potential approach [TT] to the case of a 
tightly-bound dimer being scattered off a delta poten- 
tial is particularly interesting: contrary to higher particle 
numbers A^ > 3 analytical calculations of the effective 
potential, even in higher orders, still are possible with 
reasonable effort and can be compared with numerical 
results. These can be obtained via discretisation of the 
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Schrodinger equation. The scattering of bright solitons of 
delta function like potentials has also been investigated 
on the mean-field level and has been found to successfully 
describe many aspects of the dynamics [26 28 . The non- 
linear Gross-Pitaevskii equation is not investigated here 
as it does not allow superposition states which are in the 
focus of this work. 

Considering the nonlocal quantum superposition ([!]) 
in an experiment, one should always find all particles on 
the same side of the scattering potential. As by itself 
this is not yet a convincing signature of a quantum su- 
perposition rather than a statistical mixture, we will in- 
vestigate the subsequent recombination and interference 
of the two coherent parts of the quantum superposition. 
The combined measurement of both signatures would be 
a demonstration for the existence of a Schrodinger-cat 
state. Interference experiments with matter waves can 
be of two very different kinds [6]: while the first one in- 
volves the interference of two independent condensates 
[2"9"] . the other one is based on the splitting of a single 
system into two components. In the interference of two 
independent condensates |29j an interference pattern is 
visible in a single run of the experiment but will eventu- 
ally wash out by averaging over many realisations due to 
the random phase between the condensates. Here we in- 
vestigate a complementary kind of interference where the 
interference pattern will build up in a sequence of runs: 
The density distribution for the CoM X = jj YljLi x i °f 
the system is given by 



PCc,m{x) = (S (x - X)), 



(2) 



where x denotes the spatial coordinate, whereas the 
single-particle density reads 



(3) 



The paper is organized as follows: In the next section 
we will introduce the Lieb-Liniger model j34j [35] for the 
description of a one- dimensional attractive Bose gas and 
review the effective potential approach from [11] . The 
third section covers numerical investigations of the scat- 
tering of a tightly bound dimer at a delta potential and 
the resulting interference patterns as well as the appli- 
cation of the effective potential to this two-particle case. 
In the fourth section the considerations about interfer- 
ence patterns in the CoM and single-particle density are 
extended to the experimentally relevant TV-particle case. 



II. MODELS 

A. Lieb-Liniger model of a one-dimensional 
attractive Bose gas 

One-dimensional Bose gases can be described in the 
exactly solvable Lieb-Liniger model [3H EH] where the 
many-body Hamiltonian for N bosons of mass m inter- 
acting via contact interaction with coupling constant gu 
is given by 



Hn = 



N 

y- 



.2 

2m 



Xi) 



(4) 



The resulting Schrodinger equation is separable in CoM 
and relative coordinates and can be diagonalised with the 
Bethe ansatz. For attractive interaction with < 
bound states with negative energies are expected. The 
ground state of the system is the quantum soliton 



ipo(x) = C N exp\-P y 



xA 



(5) 



i<j<e<N 



with the position Xj of the particles. While there are 
set-ups for which a single experiment yields an inter- 
ference pattern 29J, in general quantum mechanical ex- 
pectation values like equations ^ and Q only make 
statements about an average of many runs of the experi- 
ment. For the CoM density ^2§, a single experiment will 
only give a point (the CoM of all atoms measured). As 
for the double-slit experiment with single photons [50] . 
the interference pattern only builds up after repeating 
the experiment many times. The same experiments can 
also be used to determine the single-particle density (3): 
here the positions of all atoms from all experiments are 
added. Here we assume adequate experimental stabil- 
ity to garantuee that the shift of the interference pattern 
from run to run is smaller than the distance of neigh- 
boring interference maxima to avoid washing out of the 
pattern: an experiment will have to combine the low tem- 
peratures of Ref. [31] (to prepare the initial state, cf. [11]) 
with the vacuum of Ref. [32] (to suppress decoherence in- 
troduced by free atoms) and the particle-number control 
available in Experiments like [33] . 



where j3 = —mgid/2h 2 
(N-iy. - ^' 2 

N 



> 



and C 



N — 



'(2/3) 



N-l 



17] . For N ^> 1 the single-particle 

density ^ for the CoM being localized around x = is 
well approximated by the mean-field density profile [36] 



Pn (x) = J dxi 



. dx N \M*)\ 2 t ( Xi + n +Xn ) 5 (xi - x) 



ill 2 



^N 2 [cosh(^iVx)]- 



(6) 



Thus, x represents the distance from the CoM of the 
soliton. Note that the wave-function ^ already is sym- 
metric with respect to any permutation of particles. The 
ground state ([5| of the Hamiltonian Q has energy 



E (N) 



1 mg 2 d 



N (N 2 



24 h 2 

and is separated by an energy gap 

|A*I =E (N-l)-E Q (N) 

_ ,ng 2 ld N(N-l) 
— 8h 2 



(7) 



(8) 
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from a continuum of solitonic fragments. The corre- 
sponding iV-particle solutions with CoM wave vector 
K = Nk, total mass M — Nm and energies 



read 



E K {N) =E (N) 



h 2 K 2 
2M 



(9) 



density profile of the ground state soliton with the barrier 
potential: 

V eS (X) = J d N xS \X~J2x^j |V^fe(x)| 2 ^^(x,). 



A' 



ipN,k (x) = -00 (x) exp | iky^ -. 



(10) 



Hence, for low CoM kinetic energies Ek{N) — Eq(N) < 
\fi\ the scattering of the BMWS at an additional scatter- 
ing potential 



V(x) 



N 



(11) 



leading to the total Hamiltonian H = Hq + V(x), will 
be elastic. In this work we assume the scattering po- 
tentials to be delta distributions V (xj) oc 5 (xj) which 
will be motivated in section IIII1 As the soliton cannot 
break into parts nonzero probability amplitudes of the 
scattered wave-function on both sides of the scattering 
potential imply the generation of a Schrodinger-cat state. 
For higher energies particle excitations are possible and 
the soliton can break up into several 'lumps' of bound 
particles. Those excited states can be found in a gener- 
alised Bethe ansatz [17]. For N — 2 they read 



cos (kx r ) — — sin (fc|x r |) 



(12) 



with CoM and relative coordinates 

K = ki + k 2 , 
h - k 2 



x 2 - Xi, k 



B. Effective potential approach 



(13) 



Adding the additional scattering potential (11) to the 
Hamiltonian Q renders the Schrodinger equation nonin- 
tegrabel via the Bethe ansatz. Nonetheless, for suitably 
smooth scattering potentials and low kinetic energies the 
CoM motion can again be solved independently of the 
relative motion in the mathematically justified effective 
potential approach [TT] . It is essentially based on the ex- 
istence of the energy gap To solve the CoM simply 
replacing the exact potential (11 1 by NV(X) is not suffi- 
cient, see section [ill B| Instead the Schrodinger equation 
for the CoM motion is given by 



H c ff 



2M 



d\ + Kfr (X) 



(14) 



where X = ^ Sj=i x j 1S the CoM of the particles and 
the effective potential is the convolution of the internal 



III. SCATTERING OF A DIMER OFF A DELTA 
POTENTIAL 

In an experiment single atoms could serve as scat- 
tering potentials, in the following mathematically sim- 
plified as delta distributions. We therefore investigate 
the scattering of a dimer off such a delta potential in a 
one-dimensional geometry and discuss signatures of the 
resulting Schrodinger-cat states. Moreover the effective 
potential approach of section |IIB| will be applied to this 
situation which is mathematically more challenging than 
the case of broader potentials. Contrary to higher par- 
ticle numbers, for two particles analytical calculations 
nonetheless still are possible and can furthermore be com- 
pared with numerical results obtained via discretising the 
Schrodinger equation. 

Within a Bosc-Hubbard Hamiltonian [371 12H] 



latti' 



= -J 



'i+i 



*£.hj (nj - 1) 



U < 0, (16) 



describing the dynamics of N bosonic particles on a lat- 
tice, it is possible to use exact eigen-states. Here, J is 
the tunnelling strength, U denotes the pair interaction, 
alj and are the usual boson creation and annihilation 

operators and fij = a^-a- is the number operator for site 
j. For a single particle on a lattice the dispersion relation 
is given by 



E c , 



-2Jcos(kb). 



(17) 



Comparing equation ( 17) with the dispersion relation for 



suitably small lattice step sizes b under consideration of 
the lower band edge — 2 J and choosing Jb 2 = ^ then 
renders a numerical investigation of the free particle mo- 
tion possible. In the case of a dimer the two-particle 
bound states [2"TH2"4"I EM |4*D] are given by 



with 



and 



'x b ^ el exp[ikb(j + e)} ,j^£ 
±exp[ikb(j + £)} ,j = l 



U 2 



X- 



16 J 2 cos 2 (kb) 



\AJ\ cos(kb) ' 



(18) 



(19) 



(20) 
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This leads to a dispersion relation 



-4J 



/ U 2 
V16J 2 



1/2 



8J 2 b 2 



{16 J 2 + U 2 ) 



1/2' 



(21) 



Analogue to the one-particle case for small values of b 
and under consideration of the offset —4 J a comparison 
with the energy 



2fl2 



E 



JV=2,soliton 



K 2 p 



2 1,2 



K 2 k 



(22) 



of a two-particle soliton according to equation ([9| then 
enables the numerics for a tightly bound dimer moving 
in a one-dimensional geometry using the dimer wave- 
functions in a tight-binding lattice. 

In the following, like in figure [l] (a) , the initial wave 
function with mean momentum fco and width a is built 
up of the two-particle eigenfunctions ( [19] ) where k should 
be restricted to < k < f3 so that the dimer cannot break 
apart for energetic reasons. We choose the upper integral 
boundary fc max to be about twice the mean momentum: 

9 (xi,x 2 ;t = 0) c</ fcmax dfcexp(- a2(fc - fco)2 ) 

x *2,fc (%i ~ x , x 2 - x ) .(23) 



A. Interference patterns 

In the low energy-regime the scattering of a dimer at 
a delta potential 



V(x 1 ,x 2 ) = v Q [6(xi) + S(x 2 )} 



(24) 



located in the middle of the lattice, can lead to the cre- 
ation of nonlocal mesoscopic superposition states where 
both particles always should be found on the same side of 
the scattering potential. To probe the existence of such 
a Schrodinger-cat state we thus use the criterion [2T] 



dxi / dx 2 \^ (xi,x 2 ;t)\* 



(25) 



and ensure that it takes on sufficiently low values on the 
order of 10~ 8 . To investigate the interference patterns 
the contrast 



C = 



(26) 



will be calculated on a suitable interval in the middle of 
the pattern. 

Figure [l] shows numerical results for the CoM density 
distribution ([2]): The initial wave function in figure [I] 
(a) is implemented according to ( J23J ). Scattering at the 
delta potential then leads to the creation of a nonlocal 
superposition, displayed in figure [I] (b). To allow for a 
subsequent recombination of the wave packets reflection 
potentials are added at the edges of the lattice. After 
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FIG. 1. CoM density (thin line) for the scattering of a dimer 
at a delta potential V(xi, x%) = vq [5(%i) + 6(12)] . The delta 
potential situated in the middle of the lattice as well as the 
reflection potentials at the edges are displayed schematically 
(thick lines). Lattice size: 801x801 lattice points, (a) Initial 
wave function with width a/3 = 18 and mean wave vector 
ko/j3 — 0.75 centred around xo/3 — —40. (b) Schrodinger- 
cat state with p m i a0 u = 0.25 • 10~ 8 and 50%-50% splitting, 
shown at th/3 2 jm = 136. (c) Reflection of the two parts of the 
quantum superposition at reflection potentials on the edges. 
Scattering potential is switched off. (d) Recombination of 
the two parts of the wave function leading to an interference 
pattern in the CoM coordinate with maximal contrast C = 1 
at thf3 2 /m = 255. The contrast is calculated using equation 
(26b on the interval x/3 G [-10, 10]. 



being reflected, cf. figure [T] (c) , an interference pattern 
with maximal contrast C — 1 arises in the CoM density 
which can be seen in figure [l] (d) . The contrast has been 
calculated on the interval xp £ [—10, 10]. 

The results for the single-particle density ^ which are 
displayed in figure [2] differ clearly from this. The contrast 
of the interference pattern is reduced to about 0.75, again 
calculated for x(3 £ [—10, 10]. 

Summarising, the scattering of a dimer off a delta 
potential can lead to the creation of a Schrodinger-cat 
state. While the subsequent recombination of the 
two coherent parts of the wave-function leads to an 
interference pattern with maximum contrast in the 
CoM density the visibility of the pattern in the single- 
particle density is reduced. These numerical results 
will in section IV be compared with analytical con- 



siderations which can be extended to the A-particle case. 



B. Effective potential approach for the scattering 
of a dimer at a delta potential 

The existence of the energy gap {[sj) allowed the cre- 
ation of the mesoscopic superposition state in III A and 
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FIG. 2. Single-particle density (thin line) for the scattering 
of a dimer at a delta potential (thick line). Same parame- 
ters and set-up as in figure [I] (a) Initial wave function, (b) 
Schrodinger-cat state with p m iaou = 0.25 • 10~ 8 and 50%-50% 
splitting, shown at at th/3 2 jm — 136. (c) Reflection of the two 
parts of the quantum superposition at reflection potentials on 
the edges (thick line). Scattering potential is switched off. (d) 
Recombination of the two parts of the wave function leading 
to an interference pattern with reduced contrast C = 0.750 
in the single-particle density, shown at th/3 2 /m — 255 and 
computed as in figure [T] 



renders the effective potential approach from section II B 
possible. For the scattering of a dimer off a delta po- 
tential it is possible to compare the effective potential 
approach even in higher order perturbation theory with 
numerical results using again the tight-binding wavefunc- 
tions of section IIII Al 

The exact potential ( p4"| can be expressed as 

(27) 



V(xi,x 2 ) =Vx(X) + AV(x u x 2 ) 
where V X (X) = 2v Q S (A) and 

AV{x x ,x 2 ) = v [-25 (A) + S ( Xl ) + 5 (x 2 )} 



(28) 



This perturbation AV(xi,x 2 ) to the CoM contribution 
Vx(X) induces an energy shift AE on the CoM motion. 
The Schrodinger equation for the CoM wave-function 
*&(X) then is given by 



mm*)) = 



2M 



d% + V x (X) + AE(X) 



MX))- 
(29) 

The energy shift AE on the CoM motion due to the 
relative motion can be computed in Rayleigh-Schrddinger 
perturbation theory. In first order it is given by 



AEU(X) = {^ =2 \AV(x 1 ,x 2 )W 
with the ground state soliton (JsJ for two particles, 

' i0 \ =2 ( Xl ) = y/2peM-P\xA), (31) 



>(o) 

|,JV= 



(30) 



0,N=2\ 



using CoM and relative coordinates (13 1. As shown in 



|A| the energy correction adopts the form 

AE (1) (X) = -2v 5(X) + 4/3w exp(-4/3|A|) (32) 
such that up to first order perturbation theory the ef- 



fective potential in the Schrodinger equation ( 14 1 for the 
CoM can be identified as the convolution 



V^(X) + V C ^(X) 



(1)/ 



V x (X) + AEW(X) 
= 4f3v exp(-4 ( 8|X|) 



(33) 



of the scattering potential (24) with the density profile 



(31 1 of the dimer. Taking into account excited states [T7] 
gives the higher order corrections 



V$(X) 



2vim 



and 



^ exp(-4/3|X|) [1 - 80\X\ exp(-4/3|X| ;j 

(34) 



V$(X) = 



16i; 3 f r exp(-4/3|X|) 
x [C 1 (X) + C 2 (X)exp(-il3\X\) 
+C 3 (X)exp(-8^|X|)] 



(35) 



A) 



4rf 



where Ci(A) = 
-(5(A)) and 



l 

16/3 2 



to the effective potential (cf. 
T ^ + ^S(X),C 2 (X) = -(\X\ 

C 3 (X)= (8|A| 2 /3+|A| 

In the following the scattering of a dimer at a delta 
potential will be investigated numerically for suitable pa- 
rameters, on the one hand using the exact potential (24) 



8fl 



and on the other hand the results for effective potential. 
The scattering potential is again located in the middle 
of a lattice with spatial size L where —L/2 < X < L/2 
and the initial wave-function is chosen according to ( |23[ ) 
to ensure the dimer not to break apart. Figure [3] shows 
the reflection coefficient 



R(vo) 



L/2 



dX|*(X)| 2 , 



(36) 



which is increasing with the strength of the scattering 
potential vq. Neglecting the correction AE{X) leads to 
full transmission over the whole range of values for the 
height of the scattering potential and thus in zeroth order 
the dynamics is not described correctly at all. As can be 
seen in figure [3] the reflection coefficient R$(vo) for the 
exact potential is well-approximated in third order per- 
turbation theory. In the inset it is shown that already in 
first order the dynamics are qualitatively well-described. 
Especially in the experimentally interesting regime with 
R w 0.5 a convergence towards Rs(v ) can be observed: 
the second order leads to a clear improvement while the 
additional contribution in third order is small. For larger 
potential heights with R > 0.8 taking into account only 
up to 2nd rather than 3rd order contributions leads to 
better results. 
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FIG. 3. Reflection coefficient R(vo) for the scattering of 
a dimer at a delta potential V{x%,X2) = vq [S(xi) + S(x2)], 
calculated on a lattice with 801x801 lattice points and — 160 < 
xj3 < 160. The initial wave function with width a/3 = 28 
and mean momentum feo//3 = 0.375 was located at xo/3 = 
—80. Numerics for exact potential (solid line) and effective 
potential up to 3rd order perturbation theory (dash-dotted 
line). Taking only the zeroth order corrections into account 
leads to full transmission over the whole range of parameters 
and hence does not describe the dynamics correctly at all. 
Inset: Deviation AR — Rs{vo) — R E s{vo) including up to 1st 
order corrections (dotted line), 2nd order corrections (dashed 
line) and 3rd order corrections (dash-dotted line). 



IV. INTERFERENCE PATTERNS AS 
SIGNATURES OF QUANTUM SUPERPOSITION 
STATES 



In section IIII Al the interference of the two coherent 
parts of the quantum superposition has been investi- 
gated. In the following we will analytically explain the 
different values of the contrast in the CoM and the single- 
particle density and extend the two-particle results to 
considerations for N particles. We will omit the en- 
velopes of the wavepackets and restrict to the interference 
of two plane waves with opposed CoM wave vector K. If 
already in this case no interference is observed neither 
with wave-packets built up of plane waves this would be 
the case. 

For the description of a dimer we use again the relative 
and CoM coordinates (13). The wave function (10) then 
reads 



-P\x t \ p ±iKX 



(37) 



The interference of the two parts \E'+ and of the su- 
perposition 



2^e- mx ' l cos(KX) (38) 
results in the following interference pattern in the CoM- 



density distribution ([2| 

(5 (x X)) gr = f™ dx r /+~ dX|* + + \H (x X) 

oc cos 2 (xK) . (39) 

Hence, in the CoM-coordinate an interference pattern 
with fringe spacing ir/K can always be observed and the 
contrast (26) takes on its maximal value C gri coM = 1- 

For the single-particle density (J3j) the corresponding 
expectation value is 

Pgr, one — 2 ^2j = l,2^ i X — X j)) 

« xwrm t 4/?2 cos ( 2Xx ) + 4/32 + K ^ ( 4 °) 



and as the ratio f3/\K\ is related to the ratio of binding 
to CoM kinetic energy via 2/3/ \K\ - 
now is given by 



the contrast 



a 



gr,onc 



1 



ImI 



(41) 



For low CoM kinetic energies the contrast is maximal 
while for high values of -Ekin the interference vanishes: 



lim C gr , 



1. 



lim C gr! ono = 0. 



(42) 
(43) 



This agrees qualitatively with our numerical observa- 
tions, cf. figure [2j The dependence of the contrast on 
Skin /| A 4 1 is shown in figure [4] (a) . 

For CoM kinetic energies larger than the binding en- 
ergy in an experiment excitations could occur. What 
would be the influence on the visibility of the interference 
pattern in the CoM density? Repeating the same con- 
siderations as for the ground state for the excited states 
(12) leads again to an interference pattern 



lim (S (x — X )) cxc oc cos 2 (Kx) 



(44) 



with maximal contrast C CX c,CoM = 1, cf. [B] For the 
ground state we discovered the same result and looking 
separately at ground or excited wave functions additional 
phase factors are irrelevant for the contrast. However, in 
an experiment involving excitations phase factors could 
become important. Hence a superposition of the patterns 
cos 2 (Kz) and cos 2 (Kz + cf)) could lead to extinction. 

In the following we will investigate the experimentally 
relevant TV-particle case. In the low-energy regime the 
scattering of a soliton off a scattering potential has been 
predicted to lead to the creation of a nonlocal quantum 
superposition of the form ([TJ). Again the superposition 
of wave-functions ( 10 1 with opposing K leads to an in- 
terference pattern 



(Slx- 



X\ 



■ x N 



N 



)n,CoM oc cos 2 (Kx) (45) 
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with contrast Cat.coM = 1 in the CoM density distribu- 
tion while the contrast of the single-particle density is 
given by (TBI 



C 



ny/N-1 



N.one 



Skin 



sinh ( n^/N - 1 



(46) 



As shown in figure [4] (b) the contrast only takes on values 
Civ.onc ~ 1 for extremly small ratios -Ekin/l/^l and van- 
ishes otherwise. A comparison of CoM and single-particle 
density interference patterns is shown for experimentally 
realistic parameters [TT] in figure|4](c). Though the calcu- 
lations have been performed for plane waves of the CoM 
rather than wave-packets it can be deduced that also with 
wavepackets no interference in the single-particle density 
could be observed. 

Within the effective potential approach [TT] it is pos- 
sible to numerically investigate the time evolution of the 
CoM density for the scattering of a soliton consisting of 
N = 100 atoms at a delta potential which is shown in 
figure [5] The effective potential ( 15 ) being the convolu- 
tion of the scattering potential with the mean-field soli- 
ton density (|6|) is again of the form V eS ~ 1/ cosh 2 (X/2£) 
where £ = n*\gid\N is the soliton size. Clearly the split- 
ting of the wave packet due to scattering at the barrier 
potential can be observed. After switching off the barrier 
potential the coherent parts of the wave packet recombinc 
and display an interference pattern with high contrast. 

Summarising these results we have shown that the in- 
terference pattern of the CoM density distribution for 
the ground state soliton always has maximum contrast 
both for the dimer and the TV-particle BMWS. For the 
dimer the contrast in the single-particle density inter- 
ference pattern is, depending on the ratio of binding to 
kinetic energy, reduced but still visible. As shown in 
figure [4] (b) for N particles the interference vanishes for 
all but experimentally unrealistic values of -Ekin/lMh f° r 
low kinetic energies the transit time ttrans of the nonlo- 
cal superposition states in the interferometer would be 
too long in comparison to the time-scales of the main 
source of decoherence: particle losses. For the parame- 
ters of [TT] where the scattering of a soliton consisting of 
N = 100 particles is investigated a typical value would be 
^trans ^5 200 his. While one-body losses due to collisions 
with the background gas could be suppressed by a very 
high vacuum [35] three body-losses can be controlled by 
the chosen experimental parameters, leading in this case 
to a negligible loss event probability of 3%. Addition- 
ally, the very low temperatures required for the proposed 
protocol are experimentally accessible [3TJ. Besides it is 
important to note that the number of particles consti- 
tuting the soliton is experimentally well controllable: a 
number postselection of N = 100 ± 5 atoms is feasible 
[42] . So far, we have assumed the CoM of the soliton to 
be known with arbitrary precision. Nonetheless, as the 
CoM of the soliton can be determined an order of mag- 
nitude better than the soliton width [T3] the proposed 
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FIG. 4. Contrast |26| for the interference pattern in the 
single-particle density of the coherent parts of a Schrodinger- 
cat state after the scattering of a tightly bound dimer respec- 
tively a bright soliton at a delta potential, (a) Contrast | |41[ ) 
depending on the ratio of kinetic and binding energy for the 
dimer. (b) Contrast ( 46 1 for a soliton with N = 50 (solid line) , 
N = 100 (dashed lineTand TV = 200 (dash-dotted line) parti- 
cles, (c) CoM interference (solid line) of two counterpropagat- 
ing plane waves ~ exp(iKX) and *(X) ~ exp(-iKX) 

for N = 100 7 Li atoms with HK/M ~ 0.37 mm/s. In the cor- 
responding single particle density (dash-dotted thick line) no 
interference pattern can be observed for these typical experi- 
mental parameters. 



measurement of the CoM density distribution should be 
experimentally feasible. Not singling out the CoM but 
averaging over the whole measured density distributions 
would smear out this interference pattern as can be seen 
in figure [4] (b). 



V. CONCLUSION 

To summarise, for sufficiently low kinetic energies the 
scattering of tightly bound dimers or BMWS can lead to 
the creation of nonlocal mesoscopic quantum superposi- 
tion states ([I]). For the scattering of a dimer off a delta 
potential the effective potential approach from [TT] could 
be tested numerically. The approach has been found to 
describe qualitatively correct the dynamics though the 
consideration of higher order terms in the perturbation 
theory leads to clearly improved results. While so far 
only first order corrections to the effective potential have 
been investigated [TT], here we have presented a case 
where higher-order corrections become necessary. 

Furthermore, for the two-particle case we have numer- 
ically shown that the subsequent interference of the two 
coherent parts of the wave-function leads to a perfect 
contrast in the CoM coordinate while the visibility of 
the interference pattern in the single-particle density is 
reduced. Our analytical calculations confirm this be- 




FIG. 5. Time evolution for the scattering of a soliton 
consisting of N = 100 particles at a delta potential in the 
effective potential approach: numerically, by discretization 
of CoM Schrodinger equation, calculated CoM density dis- 
tribution versus dimensionless CoM coordinate X/£ and di- 
mensionless time th/£ 2 M. The initial wave packet with 
9(X,t = 0) oc exp(LKoX)exp(-X 2 /a 2 ), where a/£ = 10/3, 
is launched with an initial momentum Ko^ — 8 and is cen- 
tered around X/£ = -10. At about th/£ 2 M = 1.5 the wave 
packet reaches the scattering potential which leads to the cre- 
ation of a Schrodinger-cat state with 50% - 50 % splitting. At 
th/^M = 4 the scattering potential is switched off and the 
coherent parts recombine and interfere at th/£ 2 M « 7. Num- 
ber of lattice points: 1601. The additional reflection poten- 
tials from figures [l]and[2] are here implemented via reflection 
on the edges of the lattice. 
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Appendix A: Effective potential approach for the 
scattering of a dimer at a delta potential 



The energy change ( 30 1 in first order perturbation the- 
ory is 



AE^(X) 



-2v S(X) + 4(3v exp(-4/3|X|), (Al) 



leading to the effective potential (33). For compactness 
we omit the argument in Vx{X). 



For higher order corrections the excited states ( 13 1 for a 
dimer are required. Their relative part reads 



cos (kx r 



^sin(fc|x r |) 



(A2) 



Here the particles have been enclosed in a fictitious box 
of size L, later assumed to be infinitely large. Using 
periodic boundary conditions this leads to the norm 



l 



2L 1 



k 2 



(A3) 



In Rayleigh-Schrodinger perturbation theory with con- 
tinous eigenstates (A2) the energy correction in second 
order is given by 



haviour. For a tightly bound dimer the contrast will 
only be reduced whereas for experimentally realistic atom 
numbers on the order of N = 100 the contrast in the 
single-particle density vanishes in all but experimentally 
unrealistic regimes. Nevertheless, in the CoM density an 
interference pattern still is clearly displayed. 

Concluding, we have shown the occurrence of an in- 
terference pattern with high contrast both for the dimer 
and a soliton consisting of about 100 atoms in the CoM 
density distribution. Additionally, after the scattering 
all the particles should be clustered in a single lump, in 
repeated measurements randomly distributed on either 
side of the scattering potential. Combining these two 
characteristics can serve as a signature of nonlocal quan- 
tum superposition states. 

Note added in proof: The effective potential 
method [11] was also independently developed in [44] . 
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AE^(X) 



dkg(k) 



4 0) (2)-£f (2) 



(A4) 



where g{k) — and the argument x T in (A2) is again 
omitted. The integral is restricted to the interval < 



k < oc as the relative part of the wave-function ( A2 ) 



is symmetric with respect to k. With Eq°\n = 2) = 
-^,E^{N = 2) = >^, 



(0) 



=32/3u£yVexp(-4 / 3|X|) 



(2k\X\) - ^-sin(2fc|X|) 
k 

(A5) 



and using the Gradsteyn integrals (3.264), (3.773) and 
(3.264) we get 



AE^(X) = 



2vim 



exp(-4/?|X|) (1 - 8/3|X| exp(-4/3|X|)) . 

(A6) 

,(2) 



corresponding to the contribution \Q S to the effective 
potential in second order, given by equation (34). The 
energy correction in third order is given by 



9 



OO OO 



A ° W ^ l l [£„( 2 )-^ l( 2)][ B „(2)-^ 2 (2)] 

r ^ k* (0) IWI* (0) )l 2 (* (0) IWI* C0) ) 

£ r°° j; 1 \ 0,N = 2 I ^ x I * fcj ,rol/ I \ *0,JV = 2l Kx I *0,rel' 

~ 2^ Jo dK l 



[£ (2)-B fcl (2)] 



r 



With equations (Al), (A5| and 



4u 



cos(2fci|X|) - ^-sin(2fci|X|) 

«1 



cos(2fc 2 |A|) - ^-sin(2fc 2 |X|) 

«2 



(A7) 



as well as the Gradsteyn integrals (3.773) and (3.251) 
this gives 



E^(X) =16« 3 ^exp(-4/3|X|) 

1 



/) 4 
1 



16/3 



X 



4/3 



|X| 2 /3+|A 
1 1 



1 

8/3 



exp(-4/3|X|) 



exp(-80\X\) 



8/3 2 16/3 2 



exp(-4/3|A|) U(X) 



(A8) 



( 3) 

This corresponds to the third order correction V off to the 
effective potential, cf. (35). 



Appendix B: Interference patterns 



1. Calculation of the interference pattern for the 
excited states (12) in the two-particle case 



For the wave function ty + with K and k the corre- 
sponding wave function has been obtained by the 
replacements K — > —K and k — > —k. With the superpo- 
sition state 



cos (kx T ) — — sin (felxJ) 
k 



cos [KX) 
(Bl) 



this leads to a CoM density distribution 
lim (S (x - A)) oxc = hm / dx r / dA16AA, 



cos (fc|x r |) — — sin (fc|a; r l) 
k 



= 4 cos 2 [Kx) . 
1. iV-particle case 



cos 2 (KX) 
(B2) 



The interference pattern in the CoM density under con- 
sideration of equation is 



gij hxjy 

N 



))jV,Co 



M 



4/dxi---dxAr|^ n ,roi| 2 cos 2 (K^ N 
x S [x 



+ ---+X N ■ 



giH yx N 

N 



4 JZo d 2/ cos2 i Kx ) I dx i ■ ■■ dx N \^o, 



rcl 



= ^N 2 cos 2 (Kx) fZ, dy [cosh (SffiN (y - x))] ^ 

= AN cos 2 (Kx) . (B3) 

The interference pattern in the single-particle density ^ 
is, again using equation §6§ and Gradsteyn (3.982), given 
by 

(5 (x - xi)) N.onc 



dX cos 2 (KX) J dx x ■ ■ ■ dx N \^ 



rcl 



x 5 [X - X i+~ N +XN ) 5 (x - xx) 

cos 2 (KX) 



- 2/3 N 2 J ^ cos £°(p N ( X -x)) 

- RN 2 f°° rl y 1+cos(2K(X+x)) 
-P N J-oo dA cos ^(f}NX) 

T&T co S (2-RTx) 



2N- 



2if7T 



7r sinh (Kx) ; (B4) 

Thus, the intensity can take on the following maximal 
and minimal values 

/ max =(S(0-x 1 ))=2N+^^ w - y (B5) 

I 1 nta = (S(^-X 1 ))=2N- 

This leads to the contrast 



2Ki\ 



sinh(^)- 



(B6) 



C 



AT, one 



/9JV 



sinh 



($0 



(B7) 



Using equation ([8]) and ^ gives expression (46) for the 
contrast in dependence of 



E 



kin 



K 2 



\fi\ P 2 N 2 (N-1)' 



(B8) 
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